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We study one-dimensional reversible and number-conserving cellular automata (RNCCA) that have 
both properties of reversibility and number-conservation. In the case of 2-neighbor RNCCA, Garcfa- 
Ramos proved that every RNCCA shows trivial behavior in the sense that all the signals in the 
RNCCA do not interact each other. However, if we increase the neighborhood size, we can find 
many complex RNCCAs. Here, we show that for any one-dimensional 2-neighbor reversible parti- 
tioned CA (RPCA) with s states, we can construct a 4-neighbor RNCCA with 4? states that simulates 
the former. Since it is known that there is a computationally universal 24-state 2-neighbor RPCA, we 
obtain a universal 96-state 4-neighbor RNCCA. 

1 Introduction 

A reversible cellular automaton (RCA), and a number-conserving cellular automaton (NCCA) are kinds 
of abstract spatiotemporal models that reflect physical properties of reversibility, and conservation (of 
mass, energy, etc.), respectively. A reversible and number-conserving cellular automaton (RNCCA) is 
thus a model that has both these properties. Though an RNCCA is a very restricted subclass of a CA, its 
behavior can be complex if we increase the neighborhood size and the number of states. 

So far, NCCAs have been extensively studied, and various properties and characterizations of them 
have been given H][3l|4l|5l|7l|9l. In iflOl Moreira investigated universality and decidability of NCCAs. As 
for RNCCAs, Schranko and de Oliveira lfl31 made an experimental study on one-dimensional RNCCAs, 
and showed that an s-state n-neighbor RNCCA rule can be decomposed into s-state 2-neighbor RNCCA 
rules when s and n are very small. Garcfa-Ramos [6] proved that, in the 2-neighbor case (i.e., radius 
1/2), every RNCCA is a shift-identity product cellular automaton (SIPCA). An SIPCA is an RNCCA 
composed of "shift CAs" in which all signals are right-moving, and "identity CAs" in which all signals 
are stationary. Hence, in general, an SIPCA has both right-moving and stationary signals, but each signal 
is independent to others. Namely, every right-moving signal simply goes through stationary signals 
without affecting them. Therefore, all the 2-neighbor RNCCAs show trivial behaviors in the sense that 
the signals do not interact each other. On the other hand, Imai, Martin and Saito [ 8 ] showed that, in the 3- 
neighbor case (i.e., radius 1), there are RNCCAs in which some signals can interact with others, and thus 
they show some nontrivial behavior. However, it is not known whether there exists a computationally 
universal 3-neighbor RNCCA. 

In this paper, we investigate the 4-neighbor case (i.e., radius 3/2), and prove there is a computationally 
universal RNCCA. We show that for any given 2-neighbor s-state reversible partitioned CA (RPCA) we 
can construct a 4-neighbor 4s-state RNCCA that simulates the former. On the other hand, it is known 
that there is a computationally universal 2-neighbor 24-state RPCA lTT2"Tl . which can simulate any cyclic 
tag system proposed by Cook (H. By this, we can obtain a universal 4-neighbor 96-state RNCCA. 
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Computational universality in a variant of a one-dimensional RNCCA was studied by Morita and 
Imai [14], but this CA was not in the standard framework of NCCAs, because a partitioned CA (PC A) 
was used as an NCCA model. Though a reversible PCA is a subclass of a standard RCA, the number- 
conserving property of a PCA is somehow different from an NCCA, since each cell of a PCA has several 
parts. Namely, while each cell of a usual NCCA has a single number, that of a number-conserving PCA 
has a "tuple" of numbers. Therefore, this paper gives the first universality result of an RNCCA in the 
standard framework of an NCCA. 

2 Preliminaries 

A one-dimensional cellular automaton ( CA) is a system defined by 

A = (Z,Q,N,f,#). 

Here, Z is the set of all integers where cells are placed. Q is a non-empty finite set of states of each cell. N 
is a neighborhood, which is an element of Z m (m= 1,2, . . .). Hence N can be written as N = (n\,. . . ,n m ) 
where £ Z (i £ {1, . . . ,m}). /: Q m — > Q is a local function that determines a state transition of each cell 
depending on the states of its m neighboring cells. # € Q is a quiescent state that satisfies /(#, ...,#)=#. 
If f(qi , . . . , q m ) = q holds for qi, . . . ,q m ,q € Q, we call this relation f(qi,...,q m )=qa transition rule 
of A. Thus / can be described as a finite set of transition rules. 

A configuration over Q is a mapping a : Z — >■ Q. Let Conf(£2) denote the set of all configurations 
over Q, i.e., Conf(<2) = {a \ a : Z — > Q}. A configuration a is called finite if the set {x \x € Z Aa(x) ^ #} 
is finite. Otherwise it is infinite. The set of all finite configurations is denoted by Conffi n (<2). Applying 
the local function / to all the cells in Z simultaneously, we can obtain a global function F of A that 
determines how a configuration changes to another. More precisely, F : Conf(<2) — > Conf(<2) is defined 
by the following formula. 

Va eConf(2), ieZ: F(a)(x) = f(a(x + n\), . . . ,a(x + n m )) 

It means that the next state of a cell at the position x is determined by the present states of m cells at the 
positions x + n\, . .. ,x+n m using the local function /. If N = (— r, — r+ 1, . . . ,0, . . . ,r — l,r) for some 
natural number r G N, then A is called a CA of radius r. If N = (— r, — r + 1, . . . , 0, . . . , r — 2,r — 1) for 
some positive integer r G N — {0}, then A is called a CA of radius (2r — l)/2. Hereafter, we call a CA of 
radius r ((2r— l)/2, respectively) by a (2r+ l)-neighbor CA (2r-neighbor CA). 

Let A be a CA, and F be its global function. A is called a reversible CA (RCA) iff F is an injection, 
i.e., it satisfies the following condition. 

Vai,a 2 G Conf(£2) : ai ^ a 2 =3- F(ai) ^ -F(a2) 

A more detailed description on the definition of an RCA is found e.g. in iPTTl . 

A number-conserving CA is a one such that each cell's state is an integer, and their sum in a configu- 
ration is conserved throughout the evolving process. So far, several definitions and characterizations have 
been given for number-conserving CAs Hl[3l|4][7]]. Durand, Formenti and Roka [3 ] proved that the three 
notions, periodic-number-conserving, finite-number-conserving, and number-conserving (for infinite 
configurations), are all equivalent. In this paper, we employ the notion of finite-number-conserving to 
define a number-conserving CA. 
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Let A = (Z,Q,N, f, 0) be a CA, where Q = {0, . . . , s - 1 } (s € N - {0}), and F be its global function. 
The CA A is called finite-number-conserving, if the following condition holds. 

Va e Conf fin (£) : £ a(x) = £ F(a)(x) 

xeZ xeZ 

A CA is called a number-conserving cellular automaton (NCCA), if it is finite-number-conserving. 

A CA that satisfies both reversibility and finite-number-conserving conditions is called a reversible 
number-conserving CA (RNCCA). 

Next, we give a definition of a partitioned CA, because, in the next section, we will show a method 
of converting a reversible partitioned CA into an RNCCA. A one-dimensional partitioned cellular au- 
tomaton (PCA) is defined by 

P = (Z, (gi,. • • , Q m ), (ni, • . .,n m ),f). 

Here, Qi (i = 1 , . . . , m) is a non-empty finite set of states of the i-th part of each cell, and thus the state 
set of each cell is Q = Q\ x • • • x Q m . The m-tuple (ni , . . . ,n m ) £ Z m is a neighborhood, and f : Q ^ Q 
is a local function. 

Let pr ( : <2 — ^ <2i be the projection function such that pr,(^i, . ..,q m ) = qi for all (qi,... ,q m ) 6 2- The 
global function F : Conf(<2) — > Conf (2) of P is defined as the one that satisfies the following formula. 

Va € Conf(£2),* € Z:F(a)(x) =/(prj(a(^ + «i)),...,pr m (a(^ + ?i m ))) 

By above, one-dimensional PCA of radius 1/2 is defined as follows. 

p = (z, re,*), (o,-i),/) 

Each cell has two parts, i.e., center and right parts, and their state sets are C and R. The next state of a cell 
is determined by the present states of the center part of this cell, and the right part of the left-neighbor 
cell (not depending on the whole two parts of the two cells). Figure Q] shows its cellular space, and how 
the local function / is applied. Note that, here, the neighborhood is (0,-1) rather than (—1,0). 

It is easy to show the following lemma that states the equivalence of local and global injectivity of a 
PCA |I3 . 

Lemma 1 Let P = (Z,(Qi,...,Q m ),(n\,...,n m ),f) be a PCA, and F be its global function. Then, the 
local function f is injective, iff the global function F is injective. 

A PCA with an injective local function is thus called a reversible PCA (RPCA). 



C R C R C R 



t 


















t + i 

















x — 1 X X + 1 



Figure 1 : Cellular space of a one-dimensional 2-neighbor PCA, and its local function /. 
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3 Converting an RPCA into an RNCCA 

Lemma 2 For any given one-dimensional 2-neigbor RPCA P = (Z,(C,R), (0,-1),/), we can construct 
a one-dimensional 4-neigbor RNCCA A that simulates P and has 4\C\ ■ \R\ states. 

Proof. An RNCCA A that simulates P is given as follows. 

A = (Z,fi,(-2,-l,0,l),/,0), 

where <2 = {0,1,...,4|C| |/?|-1}. 

We need some preparations to define /. Let C,C,C,R,R, and R be as follows. 

C = {2k\R\\k = 0,\,...,\C\-\} 

C = {2(k+\C\)\R\ |Jfc = 0,l,...,|C|-l} 

c = cue 

R = {jfe|jfe = 0,l,...,|/J|-l} 

R = {k+\R\ \k = 0,l,...,\R\-l} 

R = RUR 

Each element c G C (f G R, respectively) is called a heavy (light) particle, which is a stationary (right- 
moving) particle in A as explained later. The number c (r, respectively) itself can be regarded as the 
mass of the heavy (light) particle. (Note that readers may think it strange that the particle of mass is 
both heavy and light. Though the mass could be considered as non-existence of a particle, we employ 
here the above interpretation for simplicity.) Clearly every element q G Q is uniquely decomposed into a 
heavy particle and a light particle, and thus the following holds. 

VqeQ,3c£C,3r£R(q = c + r) 

Vci,£?2 €.C,Vh,h £R (ci+n =c2 + h => c\ =c 2 /\h = f 2 ) 

We can thus regard each cell of A has a heavy particle and a light particle. We define the following 
functions pc'-Q—tC, and Pr'.Q—^R, which give a heavy particle, and a light particle associated with 
a given q 6 Q. 

VqeQ,VceC (p c (q)=c <=> 3reR(q = c + r)) 
Vq € <2,Vf G R (p R (q) = r o 3c£C{q = c + r)) 

A pair of heavy particles (c,c) G C xC (light particles (r, r) GRxR, respectively) is called a complemen- 
tary pair, if c + c = 2{2\C\ — l)\R\ (r + r = 2\R\ — 1). In the following, a complementary pair of heavy 
(light, respectively) particles is used to simulate a state in C (/?). A pair of states (<?i,<?2) £ Q 2 is called 
balanced with respect to heavy (light, respectively) particles if (pc{q\) iPcfe)) ((Pr(<1i) > Psfe))) is a 
complementary pair. The set of all balanced pairs of states (^1,^2) with respect to heavy (light, respec- 
tively) particles is denoted by Be (Br), i.e., 

B c = {(qi,q 2 ) G<2 2 I {pc{qi),Pc{qi)) € C x C 

Apc(qi)+Pc(q2)=2(2\C\-\)\R\}, 

B R = {(q u q 2 ) G Q 2 \ {p R {qi),p R (q 2 )) e R x R 

Ap R (q 1 )+p R (q 2 )=2\R\-l}. 

It is easy to see that, for any a G Conf(<2) and x G Z, the following relations hold. 



(cc(x),cc(x+ 1)) G Be => (ce(jc-l),a(jc)) 0B C A (a(jc+l),a(jc + 2)) 0B C 
(ft(jc),ft(jc+l))€fijj ^ (a(jc- l),a(jc)) 0B« A (a(x+ l),a(x + 2)) £B R 
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We now choose bijections <f>c : C — ^ C, and q>R : 7? — > R arbitrarily, and fix them hereafter. Then 
define the bijections <p c : C — > C, <p# : F — > F, <p : C x F — > £2, and <p : C x F — ► 2 as follows, where 
Q = {c + r j cGC,r GF}, and g = {c + r | cGC,f GF}. 

VceC(fc(c)=2(2|C|- -fc(c)) 
VrG*(9 s (r)=2|*|-l-^(r)) 
Vc G C, Vr G F (<p(c,r) = <p c (c) + <p R (r)) 
Vc G C, Vr G F (<p(c,r) = <p c (c) + <p R (r)) 

Now, / : Q 4 — > <2 is defined as follows. 

f Hf(^c l (Pc(qo)),9R\pR^-i)))) if (?-i,?o) e^A(go,?i) es c 

f(q-2,q-i,qo,qi)= < <p(/(<p c 1 (pc(<7-i)), <Pr '(MM))) if (<7-2,<7-i) G Br A (<7-i,<?o) G F c (1) 
I PC (<?o) + PR ) elsewhere 

Let F be the global function induced by /. For any configuration a G Conf((5), and for any y G Z, the 
value F(a)(y) is as follows. If (S(y— l),S(y)) G F s A (&(y),6;(y+l)) G Fc, tnen 

F(&)(y) = 0(/(^ 1 (pc(a(y)))^R" 1 (pR(a(y-i))))), (2) 

F(a)(y+1) = <p(/(<p c 1 (M«(3'))),<Pr 1 (pr(«(3'-1)))))- (3) 

It means the complementary pairs (/5r(&(j — 1)), pR(&(y))) and (pc(a(y)), p"c(&(y+l))) interact each 
other, and the state transition of the RPCA P is simulated. Thus, the new complementary pair of heavy 
particles (pc(F(a)(y)), pc(F(&)(y + 1))) is created at the same position as before, while the pair of 
light particles {pR{F(a)(y)),pR(F(a)(y + 1))) appears at the position shifted rightward by one cell. 
On the other hand, if ->((&{y-2), a(y- 1)) G B R A (a(y-l), a{y)) G B c )) A -.((a(y- l),a(y)) G 
B R A (a(y),a(y+l)) GF c )),then 

F(a)(y) = ^ c (a(y))+p s (a(y-l)). (4) 

The above means the light particle pR(&(y — 1)) simply moves rightward without interacting with the 
stationary heavy particle pc(ct(y)). From ©-dU), it is easy to see that the following holds for all x G Z. 

(&(*),&(*+ 1)) Gfi c <^ (F(a)(x),F(a)(x+ 1)) G Be (5) 
(&(*),&(*+ 1)) Gfifi (F(a)(jc+l),F(a)(x+2)) GS R (6) 

First, we show that A can simulate P correctly as described below. After that, we will show A is 
an RNCCA. We define a mapping f : Conf(C x R) — > Conf (Q) as follows, where a G Conf(C x R) and 

iez. 

f(a)(2x) = <p(a(x)) (7) 
f(a)(2x+l) = <p{a{x)) (8) 

The configuration a of F is thus represented by z (a) of A (see Fig. 0). We can see (pc(f(oc)(2x)), 
p c {z(a)(2x+ 1))) G F c . However, (p R (f(a)(2x- 1)), p R (f(a)(2x))) g"F R for any jc G Z. Therefore, 
by the equation (@]) we have the following. 



F(f(a))(2x) = pcMa(x)))+p R W(a(x-l))) 
F(f(a))(2x+1) = p c {<f>{a{x)))+p R {q>{a{x))) 
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By above, we can observe (p c (F(x{a))(2x)), p c {F(x (a))(2x + 1))) G B c and (Pr(F(t (a))(2x- 1)), 
Pr(F(x (a))(2*))) £ and thus the following holds by © and ©, where /? c : C x R — > C and p R : 
C x /? — >• 7? are projection functions, and F is the global function of P. 



F 2 (r(a))(2x) = (p(m c \p c {(p(a(x))))^ R HpR{(p(a(x-l)))))) 

= Q(f(pc(a(x)),p R (a(x-l)))) 

= 4>(F(a)(x)) 

= f(F(a))(2x) 

F 2 (t(«))(2x+1) = ^(m c \p c ((p(a(x))))^ R \p R ((p(a(x-l)))))) 

= <p(f{Pc(oc(x)), PR (a{x-l)))) 

= 9(F(a)(x)) 

= f(F(a))(2*+l) 



Thus, each evolution step of a configuration of P is correctly simulated by A in two steps under the 
mapping f . Its simulation process is shown in Fig. |2] 
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Figure 2: A simulation process of a 2-neighbor RPCA P by a 4-neighbor RNCCA A. The configuration 
of A at time 2t is obtained from that of P at time t by the mapping f. Here, (cj,c,-) and (r,-,r,-) (i = 1,2,3) 
are complementary pairs, and thus Cj + c; = 2(2|C| — 1)|/?| and ?,- + rj = 2\R\ — 1. 



Next, we show that A is an NCCA. From the equation (fl]), we can see mass of a particle is transferred 
within a complementary pair, or simply shifted rightward, or does not change. Therefore, it is intuitively 
obvious that A is an NCCA. But, here we show that A has the finite-number-conserving property. First, 
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from ©-(HI), we can derive the following. 

(a(x),a(x + 1)) G B c 

=> Pc(&(x)) +p c (5c(x+ 1)) = pc(F(a)(x))+p c (F(a)(x+l)) 
(a(x- l),a(x)) <£B C A (a(x),a(x + 1)) 0B C 

Pc( a W)=M^(a)W) 

(a(x),a(x+ 1)) g B s 

=>■ p s (a(x))+p i? (a(x+l)) = ^(F(a)(x+l)) + j p R (F(a)(x + 2)) 
(a(x- l),a(x)) <£B R A (a(x),a(x+ 1)) 0Br 

=> pr(cc(x)) =p R (F(a)(x+l)) 

By above, for each (= 2,3, . . .), there exist n$,n\ and «2 such that G {n,n — 1} (i = 0, 1,2), and 
«3 G {tz — l,n — 2} that satisfy the following relations. 

£ p c (a(x)) = £ Pc(^(a)W) (9) 

x=-n x=-n 
«3 "3 + 1 

£ /*(&(*)) = £ p R (F(a)(x)) (10) 

Let a G Conffi n (2). Then, the following equation holds by (© and (fTOl . 

Va G Conf(G) : £ a ( x ) = L ^( a ) (*) 

jweZ xeZ 

Therefore, A is finite-number-conserving, and thus an NCCA. 

Finally, we show A is reversible. On the contrary we assume it is not. Thus, there are two configura- 
tions a,\ , a,2 G Conf((5) such that CC\ ^ a 2 and F(ai) =P(ci2). First, we note the following. 

Vxg Z ((&i(x),ai(x+ 1)) GB C <=> (a 2 (x), a2(x+ 1)) gB c ) 
Vx G Z ((ai(x),ai(x+ 1)) G B s 44> (a 2 (i),a 2 (i+l))G5 ff ) 

If otherwise, F(cCi) ^ F(a,2) holds by the relations (f5]) and (O, and it contradicts the assumption. Since 
ai / a,2, there exists xo G Z such that pc(&i(*o)) 7^ Pc(«2(xo)) or Pi?(&i(xo)) / /5r(«2(xo))- Here, we 
prove it only for the case pc(0C\ (xo)) ^ Pc(&2( x o))> since the case Pr(&\ (x )) / /5 s (a 2 (x )) is similarly 
proved. There are three subcases: 

(i) (& i (xo-l),a i (x )) GB R A (a J -(x ),a i (xo + l)) GB c (« = 1,2), 

(ii) (a,-(x - 2), a,-(x - 1)) G B R A (a,(x - 1), a,-(x )) G B c (i = 1,2), and 

(iii) Other than the cases (i) and (ii), i.e., -i((o,(xo — l),a,(xo)) G Br A (a,(xo), S,(xo + 1)) G Be) A 
-n((a,-(x -2),a,-(x -l)) GB R A (a,(xo-l),a,-(xo)) GB C ) (i = 1,2). 

The case (i): By (f2]), the following relations hold. 

F(ai)(x ) = ^{f(q> c i {pc{oci{x ))),q> R l (p R (a l (x Q -l))))) 
F(a 2 )(x Q ) = (p(f(q> c l {p c {oc 2 {x ))),q) R 1 (p R (a2(xo-l))))) 
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From the facts pc(&i(xo)) ^ Pc(&2(xo)), <P,<Pc and <p# are bijections, and / is an injection (because P 
is a reversible PCA), F(&i)(jco) ^ F(6^)(jco) follows. This contradicts the assumption. The case (ii): 
Since it is similar to the case (i), we omit the proof. The case (iii): By (@]), the following relations hold. 

F(&i)(xo) = pc{ai(x )) + p R (ai(x - I)) 
F(a 2 )(xo) = Pc(6c2(x )) + Pr(& 2 (x - I)) 

Again F(&i)(xo) / F(&i)(xo), because pc(&i(xo)) Pc(&z(xo)), and this contradicts the assumption. 
By above, we can conclude that A is a reversible NCCA. This completes the proof. □ 

It has been shown that there is a universal one-dimensional 2-neighbor 24-state RPCA [12J. This 
RPCA can simulate Any cyclic tag system proposed by Cook [2] can be simulated by this RPCA with 
infinite but ultimately-periodic configurations. 

Proposition 1 H12II There is a computationally universal one-dimensional 2-neighbor 24-state RPCA. 
From Lemma[2]and Proposition [H the next theorem is derived. 

Theorem 1 There is a computationally universal one-dimensional 4-neighbor 96-state RNCCA. 

In lfl2l . it is shown that there is a 2-neighbor RPCA that directly simulates a given reversible Turing 
machine. Therefore, we can also construct a 4-neighbor RNCCA that directly simulates a reversible 
Turing machine. In this case, the RNCCA has ultimately periodic infinite configurations, though the 
configuration of the simulated Turing machine is finite. 



4 Concluding remarks 

In this paper, we proved that any given 2-neighbor RPCA P can be simulated by a 4-neighbor RNCCA A. 
Thus computation-universality of a 4-neighbor RNCCA is concluded in spite of the strong constraints of 
reversibility and the number-conserving property. When A simulates P, a configuration a of P is kept by 
f (a) of A as shown in Fig. [2 But, there is no need to define f as given in the equations ((TJ) and ([8]). The 
simulation works well if we use, e.g., the following f , where each two-cell-block containing <p(a(x)) 
and q>(a(x)) is separated from the next block by (k — 2) 0-state cells (k = 3,4, . . .). 

z'(a)(kx) = (p{a(x)) 
z'(a)(kx+\) = (p(a(x)) 
x'{a){kx + i) = (i = 2,3,...,jfe-l) 

Furthermore, we can see that, even if the spacing between blocks (by 0-state cells) is non-uniform, the 
simulation process goes correctly (though state transition timing of the cells are also non-uniform). 

On the other hand, it is an open problem whether a stronger result holds, i.e., whether there is a 
universal 3-neighbor (radius 1) RNCCA. It is also left for the future study to construct an intrinsically 
universal RNCCA. 

Acknowledgement. This work was supported in part by JSPS Grant-in-Aid for Scientific Research (C) 
No. 21500015 and No. 24500017. 
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